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A FAMILY OF SEQUENCES OF BINOMIAL TYPE 


WOJCIECH MLOTKOWSKI, ANNA ROMANOWICZ 


Abstract. For delta operator aD — 61)^+^ we find the corresponding polynomial se¬ 
quence of binomial type and relations with Fuss numbers. In the case D — we 
show that the corresponding Bessel-Carlitz polynomials are moments of the convolu¬ 
tion semigroup of inverse Gaussian distributions. We also find probability distributions 
vt, t > 0, for which {ynit)}, the Bessel polynomials at t, is the moment sequence. 


1. Introduction 

A sequence polynomials is said to be of binomial type (see [10]) if 

degWnit) = n and for every n > 0 and s, f G M we have 

Wk{s)Wn-kit). 

A linear operator Q of the form 

(1.2) Q = ^D + ^D‘^+ + ..., 

\ j 1! 2! 3! 

acting on the linear space R[a;] of polynomials, is called a delta operator if ci 7^ 0. Here 
D denotes the derivative operator: D1 := 0 and Df^ := n-t'^~^ for n > 1. We will write 
Q = g{D), where 

(1.3) g{x) = Y\^ + + .... 

There is one-to-one correspondence between sequences of binomial type and delta 
operators, namely if Q is a delta operator then there is unique sequence {tt’n(t)}))^o 
binomial type satisfying Qwo{t) = 0 and Qwn{t) = n ■ Wn-iit) for n > 1. These Wn are 
called basic polynomials for Q. 

A natural way of obtaining a sequence of binomial type is to start with a function / 
which is analytic in a neighborhood of 0: 

( 1 - 4 ) fix) = ^x+^x^ + ^x^ + ..., 

with oi 7^ 0. Then a binomial sequence appears in the Taylor expansion of exp {t ■ f{x)), 
namely 

(1.5) exp ((./(!)) = f; 

nl 

n=0 

2010 Mathematics Subject Classification. Primary 05A40; Secondary 60E07, 44A60. 

Key words and phrases. Sequence of binomial type, Bessel polynomials, inverse Gaussian 
distribution. 

W. M. is supported by the Polish National Science Center grant No. 2012/05/B/ST1/00626, A. R. 
is supported by the Grant DEC-2011/02/A/STl/00119 of National Centre of Science. 

1 





2 


WOJCIECH MLOTKOWSKI, ANNA ROMANOWICZ 


and the coefficients of are partial Bell polynomials of 01 , 02 ,.... More generally, 
we can merely assume that / (as well as g) is a formal power series. Then / is the 
composition inverse of g: f{g{x)) = g{f{x)) = x. 

The aim of the paper is to describe the sequences of polynomials of binomial type, 
which correspond to delta operators of the form Q = aD — Then we discuss the 

special case, D — \D‘^, studied by Carlitz [3]. We hnd the corresponding semigroup of 
probability distributions an also the family of distributions corresponding to the Bessel 
polynomials. 


2. The result 


Theorem 2.1. For a 7 ^ 0, 6 G M and for p > 1 let Q := aD — bDP~^^. Then the basic 
polynomials are as follows: Wo{t) = 1 and for n > 1 




(2.1) 

In particular Wi{t) = t/a. 


(n + j — 1)\IP 


wn{t) = V 


Proof. We have to show that if n > 1 then Qwn{t) = n-Wn-i{t). It is obvious for n = 1, 
so assume that n> 2. Then we have 


aDwnit) = 


{n + i- l)\{n - ip)V 
j! (n — jp — 




j\ {n — jp — 


and 


1^1 




j=0 


j\ {n — jp — p — l)!a”+J 


Now we substitute j' := j + 1: 


[y] 


(f) = V (n + /-2)!(n-/p + p)h^' 
(j'-l)!(n-j'p-l)!a—i+i' 

If j > 1 and jp + 1 < n then 

(n + j — l)!(n — jp) (n + J — 2)!(n — jp + p) 
j\{n-jp-l)\ {j - l)\{n - jp - 1)\ 


{n + j-2)\[{n + j - l)(n - jp) -j{n- jp + p)] 
j\{n-jp - 1 )! 

n{n + j — 2)\{n — jp — 1) n{n + j — 2)\ 
j\{n-jp-l)\ j\{n-jp-2)\ 
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and if jp + 1 = n then this difference is 0. Therefore we have 


{aD - bDP+^) Wnit) = 


t 

n I — 


71—1 




+ E7 


n{n + j — 2 )W 


j\ {n — jp — 


n-l-jp 


i=i 




= n 


{n + j — 2 )W 


j\ {n — jp — 2 )!a"'“^+'? 


-—r = n-Wn-i{t), 


which concludes the proof. 


□ 


Recall, that Fuss numbers of order p are given by and the corresponding 

generating function 


( 2 . 2 ) 


Bp{x) := 


n=0 


np + 1 \ x' 


n J np + 1 


is determined by the equation 

(2.3) Bp{x) = 1 + a; ■ Bp{xy. 

In particular 

1 — y/1 — Ax 


(2.4) 


B2{x) = 


2 x 


For more details, as well as combinatorial applications, we refer to |1]. 

Now we can exhibit the function / corresponding to the operator aD — bDP~^^ and to 
the polynomials fl2.ip . 

Corollary 2.1. For the polynomials 112.1\} we have 

exp ((. f(x)) = f; 


72=0 


n\ 


where 

(2.5) 


Proof. Since / is the inverse function for g{x) := ax — bxP~^^, it satisfies af{x) = x + 
bf{xy~^^, which is equivalent to 


af(,x) 

X 


1 + 


bxP 



p+i 


Comparing with fl2.3p . we see that af{x)/x = Bp+i (6a;^/a^+^). 

Alternatively, we could apply the formula a„ = w'.^{Qi) for the coefficients in fll.ip . □ 
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3. The operator D - 

One important source of sequences of binomial type are moments of convolution 
semigroups of probability measures on the real line. In this part we describe an example 
of such semigroup, which corresponds to the delta operator D — \D^- For more details 
we refer to ill El 0 EQ] and to the entry A001497 in [TT] . 

In view of fl2.ip . the sequence of polynomials of binomial type corresponding to the 
delta operator D — is 

+ ^ (2 n-fc-l)!t^ 

^ (n-A;)!(fc-l)!2 — 

with wo{t) = 1. They are related to the Bessel polynomials 


where Ky{z) denotes the modihed Bessel function of the second kind. Namely, for n > 1 
we have 

(3.3) Wn{t) = (1/t) = 

for n > 0. 

Applying Corollary 12.II and fl2.4p we get 


(3.4) fix) = X ■ B2{x/2) = 1 — Vl — 2x = 

n\ 

n=l 

where a„ = {2n — 3)!!. The function f{xi) admits Kolmogorov representation (see 
formula (7.15) in [12]) as 


(3.6) 


1 — \/l — 2a:i = xi + 


^H-oo 




6 

1 — wxi) , du, 

’ V2^ 


with the probability density function ^/ue on [0, +cxd). Therefore 

0 (x) = exp ^1 — \/l — 2xij 


is the characteristic function of some inhnitely divisible probability measure. It turns 
out that the corresponding convolution semigroup contained in the family of 

inverse Gaussian distributions, see 24.3 in [1]. 


Theorem 3.1. Fort > 0 define probability distribution fit := pt{u)du, where 


(3.6) 


Pi(«) := 




V2 


TTm 


for u > 0 and ptiu) = 0 for u < 0. Then {pt}t>o 'is a convolution semigroup, 
exp {t — t\/l — 2xi) is the characteristic function of pt, o,nd {'Wnit)}'^^^, defined by 
(EJP, is the moment seguence of pf 
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Proof. It is sufficient to check moments of fit- Substituting u := 2v and applying 
formula: 


(3.7) K,(t) -1(f) I exp - L 

(see (10.32.10) in [S]), we obtain 

tp r 


\ dv 


J 


„oo t.exp(^^ 
u - ^ 


'0 




du = 




exp 




2u J 


tP2 


tcyn—l fOO 


TT 


exp ^ dv = 


(t 

2 


TT 


n-1/2 




which, by fl3.3p . is equal to Wn{t). □ 

4. Probability measures corresponding to the Bessel polynomials 

In this part we are going to give some remarks concerning Bessel polynomials fl3.2p . 
First we compute the exponential generating function (cf. formula (6.2) in |5]). 

Proposition 4.1. For the exponential generating function of the seguence {yn{t)} we 
have 


(4.1) 

Proof. By 


^ yn{t) ^n ^ 


exp (j — fy/1 — 2tx) 


n=0 


ni 




we have 


E 


yn{t) 


= 


t^Wn(l/t) 


n=0 


n\ 


F^'*~^Wn+l(l/t) „ _ V UJn{i , , „ 

^ n\ ^ ^ in — 1 )! ^ 

n=l ^ ' 


n=0 


— (E 

dx ' ^ 


f^Wnil/t) 


d 


-X 


vn=0 


n\ 


dx 


= — exp (-\/l — 2tx ) , 


t t 


which leads to fl4.ip . 

Now we represent the Bessel polynomials fl3.2p as a moment sequence. 
Theorem 4.1. For n > 0 and t > 0 we have 


□ 


(4.2) 


VnW = 


"oo exp 

,71 


U 


-(77-lP 

2tu 


\/2TTtU 


du. 


Proof. Substituting u := 2tv, applying fl3.7p and 03.211 we get: 


poo exp 

.'P' 


u 


2tu 


du = 


P/t r°° 


y/2ntu ^/2nt 

gl/i( 2 ^)n+l /2 roo 


u ' 


—u 


-i/2exp( du 


\/2TTt 


= eiA 


u” 1/2 _ 


APv 


2t 2tu 
dv 



Tit 


K_n-ll2 (l/t) = yn{t), 


which concludes the proof. 

Alternatively, we could apply Theorem 13.11 and relation 03.3p . 


□ 
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Denote by vt the corresponding probability measnre, i.e. 


(4.3) 


exp 


t't := 


2tu 


V^ntu 


X(0,+oo)(m) du. 


Althongh ^ convolntion semigroup, we will see that every Pt is inhnitely 

divisible. 


Theorem 4.2. For t > 0 we have 


(4.4) = Dt/Ui/t, 

where 'yt denotes the gamma distribution with shape 1/2 and scale 2t: 

6 Xp ( ^ 

(4-5) 7 i = —^==^X(o,+oo)(m) du 

and Dtfii/t is the dilation o/pi/t by t: 


(4.6) 


Dt/ii/t 


exp 


( -P-id A 

y 2tu j 


V^ntu^ 


X(0,+oo)(m) du. 


In particular, Ut is infinitely divisible. 


Proof. From Theorem 14. II we see that the characteristic function of up 

exp (4 — j\/l — 2txi) 


(4.7) 

is the product of 


fifix) = 


\/l — 2txi 
1 


yjl — 2txi ‘ 


the characteristic function of 7 * and 

exp 


^ - 2tei ) , 


the characteristic function of which proves fl4.4p . Since both 7 ^ and are 

inhnitely divisible, so is their convolution Vf D 


Let us list some interesting integer sequences which arise from the polynomials fl3.ip 
and fl3.2p . together with their numbers in the On-line Encyclopedia of Integer Sequences 
m and the corresponding probability distribution. For their combinatorial applications 
we refer to HU: 


(1) A144301 

(2) A107104 

(3) A043301 

(4) A080893 

(5) A001515 

(6) A001517 

(7) A001518 

(8) A065919 


tc„(l), moments of pi, 

Wn{2), moments of /i 2 , 

tc„+i( 2 )/ 2 , moments of the density function u ■ p 2 {u)/ 2 , 

2 "^ • w;„(l/ 2 ), moments of the density function pi /2 {u/2) / 2 , 
|/n(l), moments of z/i, 
yn{2), moments of z/ 2 , 
ynifi), moments of z/ 3 , 

|/n(4), moments of z/4. 
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